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The two-particle decays of — > i/jir + (p + ) and — > rj c Tx + (p + ) are considered in a way taking 
into account a soft binding of quarks in the heavy quarkonia and a hard gluon exchange between 
the constituents at large recoil momenta of tp(r]c). An approximate double enhancement of the 
amplitudes is found because of the nonspectator t-channel contribution. 
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I. INTRODUCTION 

The QCD dynamics plays a significant role in an extraction of the electroweak theory parameters in the heavy 
quark sector. One of the systems allowing one to perform an exact numerical study of the heavy quark interactions, is 
the (bc) system, the heavy quarkonium with the mixed flavor. At present, the experimental search for the B c meson, 
the basic pseudoscalar lS'-state of the (be) system, is curried out at CDF Q and ALEPH j^]. 

General properties of the B c meson family can be quite reliably predicted in the theoretical investigations allowing 
one to make an objective experimental search for the B c observation (see the review on the B c physics in |||). The 
spectroscopic characteristics of (bc) family can be calculated in the framework of phenomenological nonrelativistic 
potential models 0|| and their relativistic modifications M . The strong and electromagnetic interactions conserving 
the flavor, do not give the annihilation modes of the (bc) state decays. Therefore, the excited levels radiatively 
transform into the lowest longliving pseudoscalar state decaying due to the weak interaction. The mass of this 
state, M(B+) = 6.25 ± 0.03 GeV, and its leptonic constant, fg e — 385 ± 25 MeV, can be predicted in the framework 
of potential models |HJj|, QCD sum rules and in the lattice computations |10|. The life time, t(B c ) = 0.55±0.15 



ps was estimated in several papers, where one took into account corrections caused by the quark binding inside the 
heavy quarkonium in two ways, the phenomenological one fll|l as well as in the operator product expansion for the 
weak currents of decays of the heavy quarks composing the B c meson p^ |. 

From the viewpoint of the experimental selection of the B c meson signal in a hadronic background, the preferable 
modes for the B c observation are those, wherein the final state contains the ip particle, which can be reliably identified 
in the leptonic decay, ip — > . As c-quark produced in the b — > cW* + transition, can bind the spectator c-quark of 
the B c meson with a high probability into the (cc) meson, the relative yield of ip particles in the B c decays should be 
enhanced in comparison with the branching ratio of the B u c i meson decay modes with ip in the final state. Indeed, 
under the obtained theoretical estimates in the framework of phenomenological models of the meson, one should 
expect BR(£?+ -> ipX) ~ 17 %, which is much greater than BR(B U:d -> tpX) ~ 1 %. 

As for the semileptonic decay mode of B+ — > ipl + v : estimates of its width calculated within the potential models 
fi"l| , fi"3"| and in the QCD sum rules ||, point out the essential discrepancy between results obtained in these two 
approaches (the QCD sum rule estimate of the £?+ — > ipl + v decay width is one order of magnitude less than values 
given by the different models of heavy quarkonia). As was shown in H], this deviation can be removed, if one takes 
into account the Coulomb corrections to the vertices of the meson quantum-number currents in the framework of 
QCD sum rules. 

The semileptonic mode of B c decay is suitable for the reliable experimental identification of the B c meson at a rather 
large statistics of events with B c |H[ ]. However, at the current experiments in e + e~-annihilation and hadron-hadron 
collisions, one has to expect the B c production rate, which evidently is not sufficient to identify B c in the semileptonic 
mode JlfjJlT]]. Therefore, the two-particle decay of £?+ — > ipir + allowing one to find B c practically over a single event, 
is of the greatest interest in the experimental search for B c . The estimate of its width calculated in the potential 
models, gives the branching fraction 

BR PM (B+ 7/.7T+) « 0.2% . 

However, in the semileptonic B c decay the region of low momenta for the ip particle recoil dominates, and this 
allows one to apply the approximate spin symmetry for the heavy (be) and (cc) quarkonia |L8| and reliably to use the 
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way of the transition form-factor calculation under the overlapping of the quarkonium wave functions. In contrast 
to the above transition, the two-particle modes of hadronic decays of — > ipir + (p + ) and B^ — > r) c ir + (p + ) require 
a special consideration. This is related with the fact that at large momenta of the recoil quark in the b — > cir + 
transition, the c-antiquark has to exchange by a hard gluon with the charmed c-quark being in the initial state, to 
form the bound tp(nS) or rj c (nS) state in the region of low invariant masses of the (cc) pair due to nonperturbative 
soft interactions of QCD. Thus, the feature of the two-particle hadronic B c decays is determined by the fact that in 
the 6-quark decays, the spectator quark is also heavy and, hence, at large energy release, the description of exclusive 
production of the (cc) quarkonium in the final state can not be performed in the framework of the spectator approach, 
where the quark-spectator determines only the amplitude of a soft forming of the bound state, so that the process of 
the hard weak decay can be factorized and it does not depend on the spectator. In the decays under consideration, 
this spectator picture is not valid. So, one can use the Brodsky-Lepage hard scattering formalism Jl9| . 

In this paper, we consider the exclusive B+ — > ipir + (p + ) and — > r) c ir + (p + ) decays in the framework of the hard 
approximation at large recoils with taking into account the gluon exchange to the c-quark in the initial state. In 
contrast to the spectator approach, the hard t-channel exchan ge results in the approximate double enhancement of 
the decay amplitudes, as it was recently found for _B+ — * |20| . 

In Section II we derive expressions for the amplitudes and widths of the B+ — > ipn + (p + ) and — > rj c ii + {p + ) 
decays and compare them with the spectator formulae for the b — > C7r + (p + ) transitions. Numerical estimates of the 
decay widths are given in Section III, where theoretical uncertainties of the values are discussed. The obtained results 
are summarized in the Conclusion. 



II. CALCULATION OF TWO-PARTICLE WIDTHS OF B c 



In the framework of the nonrelativistic formalism for the heavy quark binding into the S'-wave quarkonium, we 
assume that the momentum of the quark, composing the meson, is equal to Pq = itiqv^, where is the four-velocity 
of quarkonium, so that the quarks inside the meson move with the same four-velocity v. Moreover, the quark-meson 
vertex with nontruncated quark lines corresponds to the spinor matrix 

. 1 + }M nS 

Vv ~ e — w 

for the vector quarkonium with e p , being the polarization vector, and 
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v /M n c 



2^3 ' 



for the pseudoscalar quarkonium. Here M n $ is the nS-level mass and / is related with the value of configuration 
wave function at the origin 



/ = 



12 



M, 



nS 



l*n S (0)| 



The / quantity can be related with the leptonic constants of states 

{Q\M0)\V) = ifvM v e M , 
(0|J Bp (0)|P) = if PPli , 

where J^ix) and J^^x) are the vector and axial- vector currents of the constituent quarks. Then the allowance for 
the hard gluon corrections in the first order over a s |7|]2l|-|23|| results in 



/ = fv 
f = fp 
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where are the masses of quarks composing the quarkonium. For the vector currents of quarks with equal masses, 
the BLM procedure of the scale fixing in the "running" coupling constant of QCD [Q gives |23| 

of = a™(e-^m Q ) . 



2 



For the quarkonium with mi 7^ m2, we assume 



, H 



-11/12 



mim 2 ) 



Note, in the given estimates one considers the hard gluon corrections to the quark- antiquark annihilation currents. The 
corresponding factors are known exactly, and that is surprisingly, they can be obtained by the symbolic substitutions 
mi —* —mi, V «-> P from the exact expressions for the hard gluon factors of the quark-to-quark transition currents 
PH , considered in HQET pq|, at the prescription of the absolute value for the logarithm argument. However, these 
substitutions do not lead to valid evaluations of the BLM scales determining aj? . The corresponding BLM scales in 
HQET have been calculated by M. Neubert p6| ], and they do not give the exactly known result for the quark-antiquark 
annihilation vector current [j23fl . 

Further, the factor of the colour wave function <J* J '/\/3 stands in the quark-meson vertex. 

The 7r meson current corresponds to the axial- vector current of weak transition = f^p^- So, the given factor- 
ization neglects possible final state interactions, which really seem to be small (see discussion in p7y). 



7T+(p+) 



7T+(p+) 




4>M B+ 




Figure 1: Diagrams of the 5+ 
between the constituent quarks. 



V'( 7 ?c) 7r+ (p + ) decays with the hard gluon exchange 



Then the matrix elements of — > ipTT + and £>+ — *• t] c tt + decays calculated according to the diagrams in Fig. 1, 
take the form 



T(5+ - ^r+) = ^§ V bc fJ Be U 



M 2 



(e • v) ax 

2 \ ( Q J\yf 2 



T(B C _ ^ ) = _ Vbc — U B J Vc Mml(y-l) 2 01 ' 



(3) 
(4) 



where v is the four-velocity of B c meson, M is its mass, e is the polarization vector of ip particle, my, is its mass, 
y = v ■ ity is the product of the i? c and ip four- velocities, 



V 



M 2 + m 2 , 
2Mrriip 



In (||) the notations for the r\ c state are analogous to the described ones. The a\ factor is caused by the hard 
gluon corrections to the effective four-fermion weak interactions of quarks. a\ is evaluated for free quarks, so that 
<X\ = 1.22 ± 0.04 ||. In present calculations, this correction is at the same level as the hard corrections to the decay 
constants as well as the gluon propagator. The value of the QCD coupling constant is determined by the gluon 
virtuality, k 2 = — m^ (y — l)/2, and it will be discussed in next Section. 

The corresponding virtualities of c and b quarks, interacting with the hard gluon, are equal to 



;.2 



= 2k 



ml = 2kgM/m^ iVc 



Note, that numerically the virtual c-quark in the second diagram in Fig. 1 is in the ^-channel, since its four-momentum 
squared has a negative value. Therefore, one can see, that the corresponding contribution into the — > '0( 7 7c) 7r+ 
decay is definitely nonspectator, and the considered process is certainly hard. 
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Prom Eqs.(|,|) one gets the expressions for the total widths of the — > ipir + and — > ry c 7r + decays 

r(B c ) = G F |F 6c | _^/ ff / fl ^(^__j ____ ai , (5) 

, (3M 2 -2Mm fl + m 2 ) 2 
T(B+ -> r /c ^+) = r(S+ -> V^+) ■ i ^ ^ . (6) 

As for the analogous two-particle decays with p + in the final state, one uses the approximate factorization of 
the transition current of the virtual W* + boson into ir + or p + , and one finds that the only difference between the 
squares of amplitudes for the pseudoscalar and vector states of the light quark systems, is the substitution of the 
quantity f p m 2 p (—g pl , + p^p^/m 2 ) instead of the ir + meson current tensor f^p^p^. Then one can easily observe that 
after the summation over the p meson polarizations, the squares of the matrix elements coincide up to the factor, so 



T(B+ -> 4p+) _ T(B+ -> VcP +) _ I 
T(B+ - Vtt+) ~ T(B+ - Vc ir+) ~ f 2 

in the leading order over the small parameters, rn 2 /m 2 , m 2 p /m 2 . 
Write down the expressions for the spectator decays of 6-quark 



(7) 



r 6 - c P +) = r & - c-7r+) ^ 1 + /V c ' . (9) 



From Eqs.(||,||) one can see that for the spectator decays the relation between the yields of p + and ir + mesons 



r(6 - c P +) ^ /, 



T{b -> C7T+) /a 

is valid in the leading approximation over the square of ratios of the p meson mass over the heavy quark masses, as 
it takes place for the transitions between the mesons (see (0)). In the spectator decays the accuracy of the leading 
approximation used, is about 4 %, that also points out the magnitude of the correction terms to relation ([7]) for the 
mesons. 

The breaking of the spectator picture at large recoils in nonhadronic decays of B£ was also recently considered in 
pifl , where one studied the £?+ — * 7 mode due to the flavor-changing neutral current of b — ► sj. The discussion of 
the heavy quark symmetry at large recoils in the heavy- light meson transitions was given in p9| , where the peaking 
approximation was introduced. This approximation corresponds to the quark-meson vertices used in this paper. 



III. NUMERICAL ESTIMATES 



The accuracy of the given calculations is basically restricted by the uncertainty in the choice of the QCD coupling 
constant value. In Eqs.(^,^) a s can be evaluated at the scale typical for the charm quark physics a s ~ 0.30. The 
higher order corrections are beyond the scope of this paper. Nevertheless, to evaluate the possible value of these 
corrections, one can use the BLM procedure including the light quark loops in the virtual gluon propagator. So, a s 
is given at the scale of the gluon virtuality 

^ = -m^ o (y - l)/2 « -1.2 GeV 2 , (10) 

so that 

a s = a^(-e-^k 2 ) . 

As one can see from Eq.(10|), the virtuality of hard gluon is comparable with the square of charm quark mass, and 



it indicates the applicabi 



ity of the hard process factorization. Moreover, the scheme-independent value of the a s 



argument e Cecheme k 2 /A 2 cheme is quite large, and it is close to 15. Numerically, the BLM fixing of the QCD coupling 
constant gives af LM s» 0.57. In the above procedure the gluon virtualities are taken into account. However, the 
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quark virtualities can be essential, too. Indeed, the charmed quark virtuality is valued in the intermediate region, 
where it is less than the 6-quark one and greater than the gluon virtuality. So, in the following estimates we will use 
the value 

a s = off* (ml - k 2 ) pa op (-2A£) = 0.33 ± 0.06 . 

The uncertainty in the a s value appropriate for the given process indicates a possible large role of higher order 
corrections. 

In numerical estimates we suppose [ pof 

\V bc \= 0.041 ±0.003 , 
and we use the one-loop expression for the a s evolution 

-in 



a s (m ) 



where f3o( n f) — 11 — 2n//3, rif is the number of quark flavors with m nf < m, 



A( n/ ) - A (n/+1 ) 



n/ + l \ Wo(»/) 



A(„ /+ i) 



Using af s (m 2 z ) = 0.117 ± 0.005 jn), one finds that A (5) = 85 ± 25 MeV and A (3) = 140 ± 40 MeV. One estimates 

af s (ml) = 0.20 ± 0.02, that is quite reasonable. 

The fs c constant was estimated in the framework of the QCD sum rules f|-g 

f Bc = 385 ± 25 MeV, 

and it is in a good agreement with the scaling relation for the leptonic constants of 15 heavy quarkonia 

f 2 ( M \ mim 2 



const. , 



M \Hi2 J mi + m 2 

Then, the account for the hard gluon corrections gives 

fi> — frjc = 542 ± 50 MeV, (11) 
f Bc = 440 ± 40 MeV. (12) 

To calculate the branching ratios we evaluate the total B c meson width according to the formula by fjjljj 

T{B C ) w T{B) + (0.6 ± 0.1)r(D+) + T(ann.) , 

where T(B) is the contribution of 6-quark decays with the spectator c-quark, T(D + ) determines the contribution of 
c-quark decays with the spectator 6-quark and with the account for the phase space reduction, because of the c-quark 
binding inside B c (i.e. one takes into account the deviation from the exact spectator consideration), and T(ann.) is 
the contribution of annihilation channels depending on | V bc \ and fg c ■ Then 

T(B C ) = (1.2 ±0.2) • 10~ 3 eV ' 



0.55 ± 0.15 ps 



The recent estimates of the T(B C ) total width calculated within the operator product expansion approach and including 
also the bound quark effects, annihilation channels as well as the Pauli interference in the final state of B c decay, is 
in a good agreement with the given value [fl2| . 

Supposing ai = 1.22 ± 0.04 g, f p = 220 MeV, one finally finds 

F(B+ -> ^+) = (9.2 ± 2.3) • 10- 6 eV = 1 , (13) 

69 ± 17 ps 

T(B+ -> i,p+) = (24 ± 6) • 10- 6 eV = 22 ± 5ps , (14) 
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and 

BR HS (B+ -> ^tt+) = 0.77 ±0.19% , BR HS (S+ -> r; c vr+) = 1.00 ± 0.25% , (15) 
BR HS (B+ -> ^ + ) = 2.25 ±0.56% , BR HS (S+ -> ?7 c p + ) = 2.78 ± 0.70% . (16) 

Further, the purely spectator decays of b — > C7r + (p + ) have the following branching fractions with respect to the 
total B+ width 

BR spoc (6^OT+)«0.64%, (17) 
BR spcc (6^cp+)«1.8%, (18) 

The matrix element, corresponding to the first diagram in Fig. 1, is approximately equal to the matrix element, 
following from the second diagram and, hence, estimates ( pi3||l6| ) are enhanced by a factor of four due to the i-exchange 
nonspectator contribution. 

As for the n5*-excitation yields of the (cc) quarkonium in the £?+ decays, we note that the corresponding branching 
fractions are determined by the rescaling of the leptonic constants and phase spaces. So, 



f2 ]\,f2 _ m 2, 

BR HS (3+ - tf(nS>r+(,+)) = BR HS (B+ - ^ + (p + )) ^ ^ . (19) 



The experimental values of leptonic constants are in a good agreement with the scaling expression 

fnS _ 1 m i> 

/| n m nS ' 

So, neglecting the differences in the masses of ^{nS) and rj c (nS) states, one gets 

BR HS (B+ -» ip(2S)w+) _ BR HS (B+ -» r7 c (25)7r+ 



(20) 



BR HS (B^ -> ^7T + ) BR HS (S^ -> T7 C 7T+) 

for instance. The same values for the ip(nS)p + state yields can be rewritten down 



0.36 , (21) 



IV. CONCLUSION 



In this paper we have shown that in the Bj — > ipir + (p + ) and B+ — > t] c tt + (p + ) decays the large momentum of the 
recoil ip or r\ c particle leads to the fact that the formalism of the weak transition form-factor calculation, based on 
the overlapping of the nonrelativistic wave functions for the heavy quarkonia, is not valid. The hard gluon exchange 
with the spectator quark results in the large virtuality of heavy quark in the weak transition current. The amplitude 
of the weak decay with the hard exchange by gluon can be calculated in the framework of QCD perturbation theory 
and this amplitude can be factorized from the amplitude of soft binding of heavy quarks in the quarkonium. The 
calculations with the account for this hard-soft factorization result in 

BR(B+ -> %lm+) = 0.77 ± 0.19% , 
BR(E+ -> tP p +) « 2.78 ■ BR(B+ -> ^tt+) , 
BR(B+ -> VcK + (p + )) » 1.28 • BR{B+ -» ^vr + (p+)) , 

where the accuracy is basically restricted by uncertainties in the evolution scale of the "running" QCD coupling 
constant, in the c-quark mass, and the total B c width. The given estimate for the branching ratio of the B^ — > ipn + 
decay mode is significantly larger than the extrapolation results by the potential models. This value strongly enhance 
the probability of B c observation in the current Tevatron and LEP experiments with vertex detectors. 

The author thanks A.Razumov and S.Slabospitsky for the software help in the figure creation and A.A.Likhodcd 
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